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1.  Introduction 

1.1.  Parameter  identification  problems  modelled  by  PDEs 
Parameter  identification  problems 

-  arise  in  the  inverse  problems  theory 

-  are  concerned  with  the  identification  of  physical 
parameters  from  observations  of  the  evolution 

of  a  system 

f  fx,  w(m\  E(nA  =  0,  xencRp  (1) 

E  =  E(x)  is  the  parameter  and  w  =  w(x)  is  the  data 

-  assume  that  the  parameter  E  and  the  data  ware  analytical 
functions 

-  the  PDE  (4)  sometimes  augmented  with  certain  boundary 
conditions  is  called  the  inverse  problem  associated  to  a  direct 
problem.  The  direct  problem  is  the  same  equation  but  the 
unknown  function  is  the  data  w,  for  which  certain  boundary 
conditions  are  required 

•  In  general,  these  are  ill-posed  problems:  they  do  not  ful¬ 
fill  Hadamard’s  postulates  for  all  admissible  data:  a  solution 
exists,  the  solution  is  unique  and  the  solution  depends  contin¬ 
uously  on  the  given  data. 

•  Arbitrary  small  changes  in  data  may  lead  to  arbitrary  large 
changes  in  the  solution 

•  The  iterative  approach  of  studying  these  problems  is  a  functional- 
analytic  setup  with  a  special  emphasis  on  iterative  regulariza¬ 
tion  methods  [i] 


[1]  H.  W.  Engl,  M.  Hanke  and  A.  Neubauer,  Regularization  of 
Inverse  Problems,  Kluwer,  Dordrecht,  1996 


1.2  Symmetry  reductions  for  a  PDE 


Consider  a  n-th  order  PDE 

T  =0,  u  —  u(x ),  x  e  Q  c  Rp 

Symmetry  reductions 
— ►  analytical  solutions  of  the  PDE  (2) 

— ►  a  reduced  order  model  related  to  (2) 

f’(z,v^n^=  0,  V  =  v(z),  z  e  Q  C  Rp_1 

where  z  =  z(x)  and  y  =  A(x,y(x)). 

1.  classical  Lie  symmetries  —>  symmetry  group 

2.  nonclassical  (conditional)  symmetries 

3.  If  (2)  depends  on  an  arbitrary  function: 


equivalence  transformations 


1.  Classical  Lie  symmetries  -  (local)  Lie  groups  of  transforma¬ 
tions  acting  on  the  space  of  the  independent  variables  x  and 
the  space  of  the  dependent  variable  u  with  the  property  that 
they  leave  the  PDE  unchanged;  the  classical  Lie  method  ([2] 
and  [3]) 

2.  Nonciassicai  (conditional)  symmetries  -  (local)  Lie  groups 
of  transformations  acting  on  the  space  of  the  independent 
variables  x  and  the  space  of  the  dependent  variable  u  with  the 
property  that  they  might  not  leave  the  PDE  unchanged,  but 
using  them  we  can  find  exact  solutions  of  the  studied  model; 
e.g.,  the  nonciassicai  method  [4]  and  the  direct  method  [5] 

3.  Equivalence  transformations  [6]  -  classical  Lie  symmetries 
acting  on  the  space  of  the  independent  variables  x,  the  space 
of  the  dependent  variable  u  and  the  space  of  the  arbitrary 
function;  the  classical  Lie  method  applied  to  the  studied  PDE 
where  the  dependent  variable  u  and  the  arbitrary  function  are 
both  unknown  functions  [7] 


[2]  S.  Lie,  Gesammelte  Abhandlungen ,  4,  B.  G.  Teubner,  Leipzig,  1929, 
pp.  320-384 

[3]  P.  J.  Olver,  Applications  of  Lie  Groups  to  Differential  Equations ,  Grad¬ 
uate  Texts  Math.  107,  Springer  Verlag,  New  York,  1986 

[4]  G.  W.  Bluman  and  J.  D.  Cole,  The  general  similarity  solutions  of  the 
heat  equation,  J.  Math.  Mech.,  18  (1969),  pp.  1025-1042 

[5]  P.  A.  Clarkson  and  M.  Kruskal,  New  similarity  reductions  of  the  Boussi- 
nesq  equation,  J.  Math.  Phys.,  30  (1989),  pp.  2201-2213 

[6]  L.  V.  Ovsiannikov,  Group  analysis  of  differential  equations,  translated 
by  W.  F.  Ames,  Academic  Press,  New  York,  1982 

[7]  N.  H.  Ibragimov,  CRC  Handbook  of  Lie  Group  Analysis  of  Differential 
Equations,  Vol.l:  Symmetries,  Exact  Solutions  and  Conservation  Laws, 
CRC  Press,  Boca  Raton,  1994 


1.3.  Symmetry  reductions  related  to  a  parameter  identifica¬ 
tion  problem  modelled  by  a  PDE 

Consider 

T  (x,  w(m\  =  0,  X  e  n  C  Rp  (4) 

E  —  E(x)  is  the  parameter  and  w  —  w(x)  is  the  data 

-  how  parameter  identification  problems  can  be  analyzed  with 
the  tools  of  symmetry  analysis  theory  (although  the  symmetry 
analysis  theory  is  a  common  approach  for  studying  PDEs) 

1.  E  unknown,  w  given 

— >  classical/nonclassical  symmetries 

2.  E  unknown,  w  arbitrary 

— >  classical/nonclassical  symmetries 

3.  E ,  w  unknown 


equivalence  transformations 


-  how  to  incorporate  invariant  data  in  non-invariant  boundary 
conditions 

!  Q  bounded 

!  w  satisfies  certain  boundary  conditions 

•  <9Q  =  {z  =  k}  where  2  is  a  similarity  variable;  w  and  Q  share 
the  same  invariance  [8] 

Prop.  If  G  is  a  (local)  Lie  group  of  transformations  on  Rp  and 
H  :  M  — ►  R/  is  a  smooth  function,  then 

H  is  G-invariant  iff  every  level  set  {x  :  H(x)  =  k},  k  g  B!  is  a 
G-invariant  subset  of  M 

E.g.  w  =  g(z),  d£2  =  {z  =  k}  is  the  fc-level  of  z,  i.e.,  X(w)  =  0, 
X  is  the  tangent  vector  to 

— ►  rotational  symmetries 

•  dQ.  =  {z  =  k}  where  2;  is  not  a  similarity  variable;  w  and  Q 
do  not  share  the  same  invariance 

E.g.  w(x,y)  =  z(x,yYn  —  krn  and 

dQ.  =  {(x,y)  :  z(x,y )  =  k},  m  >  1  natural,  k  fixed 

-  elliptical  domains 

-  rounded  square  domains  (Lame  ovals): 
x2n  _|_  y2n  <  if  n  >  2  natural 

— ►  scaling  symmetries 


[8]  G.  W.  Bluman  and  S.  Kumei,  Symmetries  and  Differential  Equations, 
Appl.  Math.  Sci.  Vol.  81,  Springer,  new  York,  1989 


2.  A  mathematical  model  arising  in  the  car  windscreen  design 
2.1.  Introduction  to  the  manufacture  of  car  windscreens 


The  main  processes  used  in  the  manufacture  of  car 
windscreens  are  the  gravity  sag  bending  process  and 
the  pressing  process 

a)  The  gravity  sag  bending  process 

-  a  piece  of  glass  is  placed  over  a  rigid  frame  with  the  desired 
edge  curvature  and  heated  from  below 

-  the  glass  becomes  viscous  due  to  the  temperature  rise  and 
sags  under  its  own  weight 

-  the  final  shape  depends  on  the  viscosity  distribution  of  the 
glass  obtained  from  varying  the  temperature 

b)  The  pressing  process 

-  a  piece  of  glass  is  heated  and  forced  between  two  molds  or 
dies  of  the  desired  shape 

-  the  glass  is  allowed  to  cool  and  set  into  the  final  shape. 


The  two  main  considerations  during  these  processes  are  the 
cost  and  the  optical  quality  of  the  product. 

The  optical  quality  is  affected  by 

-  the  quality  of  the  glass:  the  initial  glass  used  by  manufac¬ 
turers  are  produced  to  a  high  optical  quality 

-  the  required  windscreen  shape:  the  windscreen  shapes  are 
designed  in  order  to  minimize  any  adverse  optical  effects:  for 
example,  manufacturers  try  to  avoid  shapes  with  areas  where 
Gaussian  curvature  becomes  negative,  which  can  give  regions 
of  poor  visibility 

-  defects  introduced  during  production:  it  is  important  that 
the  windscreen  does  not  contain  distortions  or  imperfections, 
so  that,  it  is  essential  that  the  production  process  to  not  affect 
the  optical  quality  of  the  windscreen 

When  a  windscreen  in  produced  using  both  of  the  two  pro¬ 
cesses,  unavoidable  defects  due  to  the  production  process  are 
introduced  into  the  glass  sheet 

This  is  mainly  due  to  stresses  caused  by  the  glass  itself 

In  order  to  minimize  these  stresses,  there  is  a  final  process 
after  the  initial  production,  the  so-called  annealing  process, 
when  the  windscreen  is  gently  heated  for  obtaining  the  desired 
shape 

Sag  bending  process  is  found  to  introduce  fewer  defects  into 
the  glass  than  pressing  process 

Windscreens  produced  by  sag  bending  require  less  annealing 
and  so  are  cheaper  and  quicker  to  produce 

!  It  is  not  obvious  which  shapes  can  be  made  using  the  gravity 
sag  bending  process 


It  has  been  shown  that  the  sag  bending  process  can  be  con¬ 
trolled  (in  a  first  approximation)  in  the  terms  of  Young’s  mod¬ 
ulus  E,  a  spatially  varying  glass  material  parameter,  and  the 
displacement  of  the  glass  w  can  be  described  by  the  thin  linear 
elastic  plate  theory  ([9],  [10]  and  [11]) 

The  model  is  based  on  the  linear  plate  equation 

( E(WXX  +  UWyy))xx  +  2(1  -  U)  (EWXy)xy 

i  o  r  i  2  \  £ 

+  (E  O yy  +  VWXX))yy  =  - h3~ -  °n  (5) 


where 


w  =  w(x,y )  is  the  displacement  of  the  glass  sheet  (the  target 
shape)  occupying  a  domain  Q  c  R2 

E  =  E(x,y )  is  Young’s  modulus:  a  positive  function  that  can 
be  influenced  by  adjusting  the  temperature  in  the  process  of 
heating  the  glass 


/  is  the  gravitational  force 


V 


is  the  glass  Poisson  ratio 


h  is  thickness  of  the  plate 


[9]  D.  Krause  and  H.  Loch,  Mathematical  Simulation  in  Glass  Technology , 
Springer  Verlag  Berlin  Heidelberg,  new  York,  2002 

[10]  P.  S.  Manservisi,  Control  and  Optimization  of  the  sag  bending  process 
in  glass  Windscreen  design,  in  L.  Arkeryd,  J.  Bergh,  P.  Brenner,  R.  Petter- 
son,  Progress  in  Industrial  Mathematics  at  ECMI98,  Teubner,  Stuttgart, 
1999,  pp.  97-105 


[11]  O.  S.  Narayanaswamy,  Stress  and  structural  relaxation  in  tempering 
glass,  J.  Am.  Ceramic  Soc.,  61  (1978),  pp.  146-152 


(E  (WXX  +  UWyy))xx  +  2(1  -  V)  ( EWXy)xy 

,  ^  12(1  —  i/2)/  ^ 

+  (£  (.Wyv  +  vwxx))yy  = - — -  (5) 


-  the  direct  problem  (or  the  forward  problem):  for  a  given 
Young  modulus  E,  find  the  displacement  w  of  a  glass  sheet 
occupying  a  domain  Q  before  the  heating  process 

-  the  PDE  (5)  is  an  elliptic  fourth  order  linear  PDE  for  the 
function  w  =  w(x,y) 

-  two  problems  related  to  (5)  have  been  studied:  the  clamped 
plate  case  and  the  simply  supported  plate  case  [12] 


-  we  study  the  clamped  case,  in  which  the  following  boundary 
conditions  are  required:  the  plate  is  placed  over  a  rigid  frame, 


w(x,y)\dn  =  0, 


and  respectively, 


dw 

dn 


which  means  the  (outward)  normal  derivative  of  w  must  be 
zero,  i.e. ,  the  plate  is  not  allowed  to  freely  rotate  around  the 
tangent  to  <9Q 


-  the  associated  inverse  problem  consists  of  finding  Young’s 
modulus  E  for  a  given  data  w  in  (5) 


[12]  E.  H.  Mansfield,  The  bending  and  stretching  of  plates ,  Cambridge 
University  Press,  1989 


(E  (WXX  +  VWyy))xx  +  2(1  -  V)  ( EWXy)xy 

,  ^  12(1 -u2)f  ^ 

+  (E  (Wyy  +  UWXX))yy  =  - — -  (5) 

-  the  PDE  (5)  is  a  linear  second  order  PDE  for  Young’s  mod¬ 
ulus  that  can  be  written  as 

(^xx  “1“  V  yy)^xx  “t-  2(1  v')  \J Xy^ xy  (  yy  “I"  V  VXx) Eyy 

+  2(Aw)xEx  +  2(Aw)yEy  +  (A2w)E=l,  (8) 
after  the  scaling  transformations  w  — >  or  E  — >  }  E,  with 

U  —  12(1— t^2)/ 

K  h 3 

-  in  (8),  A  denotes  the  Laplace  operator 

-  the  main  problem  in  the  car  windscreen  design  is  that  the 
prescribed  target  shape  w  is  frequent  s.t.  the  discriminant 

D  =  (1  -  v)2wly  -  (wxx  +  UWyy){Wyy  +  VW XX) 

of  the  equation  (8)  changes  sign  in  the  domain  Q,  so  that  (8) 
is  a  mixed  type  PDE 

-  this  is  one  of  the  reason  for  which  optical  defects  might 
occur  during  the  process 

-  the  equation  (8)  would  naturally  call  for  boundaries  condi¬ 
tions  for  E  on  <9Q  in  the  purely  elliptic  case  (when  D  <  0),  and 
Cauchy  data  on  a  suitable  (non-characteristic  part)  r  c  <9Q  in 
the  purely  hyperbolic  part  (for  D  >  0) 


There  is  a  recent  interest  in  studying  this  inverse  problem 


-  Mansfield  [12]  showed  that  a  constant  Young’s  modulus  cor¬ 
responds  to  a  data  which  satisfies  the  nonhomogeneous  bihar¬ 
monic  equation  A 2w  =  ct. 

-  Salazar  and  Westbrook  [13]  studied  the  case  when  the  data 
and  the  parameter  are  given  by  radial  functions 

-  Temple  [14]  gives  a  survey  on  this  subject  and  considered 
new  examples  of  data  for  rectangular  frames 

-  Kugler  [15],  [16]  used  a  derivative  free  iterative  regularization 
method  for  analyzing  the  problem  on  rectangular  frames 

-  Engl  and  Kugler  [17]  studied  a  simplified  model  for  the  in¬ 
verse  problem  on  circular  domains 


[12]  E.  H.  Mansfield,  The  bending  and  stretching  of  plates ,  Cambridge 
University  Press,  1989 

[13]  D.  Salazar  and  R.  Westbrook,  Inverse  problems  of  mixed  type  in  linear 
plate  theory ,  Eur.  J.  Appl.  Math.,  to  appear 

[14]  D.  Temple,  An  inverse  system:  An  analysis  arising  from  Windscreen 
manufacture,  M.  Sc.  Thesis,  Department  of  Mathematics,  Oxford  Univer¬ 
sity,  Oxford,  United  Kingdom,  2002 

[15]  P.  Kugler,  A  derivative  free  LandWeber  method  for  parameter  iden¬ 
tification  in  elliptic  partial  differential  equations  With  application  to  the 
manufacture  of  car  Windscreen  s,  Ph.D.  Thesis,  Institute  for  Industrial 
Mathematics,  Johannes  Kepler  University,  Linz,  Austria,  2003 

[16]  P.  Kugler,  A  parameter  identification  problem  of  mixed  type  related 
to  the  manufacture  of  car  windscreen  s,  SIAM  J.  Appl.  Math.,  to  appear 

[17]  H.  W.  Engl  and  P.  Kugler,  The  influence  of  the  equation  type  on 
iterative  parameter  identification  problems  Which  are  elliptic  or  hyperbolic 
in  the  parameter,  Eur.  J.  Appl.  Math.,  14,  2  (2003),  pp.  129-163 


Our  aim  is  to  find  the  symmetry  reductions  related  to  the 
PDE  (8)  hidden  by  the  nonlinearity  that  occurs  between  the 
data  and  the  parameter 

— ►  group  of  transformations  that  leave  the  equation  unchanged, 
which  also  relate  the  inverse  and  the  direct  problem 

— ►  knowledge  of  the  invariants  of  these  group  actions  allows 
us  to  write  the  target  shape  and  the  parameter  in  terms  of 
them,  and  therefore,  to  reduce  the  order  of  the  studied  model 

— ►  we  find  again  the  obvious  result  that  a  Young’s  modulus 
constant  corresponds  to  a  data  which  is  a  solution  of  an  non- 
homogeneous  biharmonic  equation 

— ►  the  circular  case  considered  by  Salazar  and  Westbrook  is  a 
particular  case  of  our  study 

— ►  other  target  shapes  which  are  not  radial  functions  can  be 
considered 

-  the  equation  (8)  is  invariant  under  scaling  transformations 

-  target  shapes  modelled  by  homogeneous  functions  can  be 
analyzed  as  well. 

-  in  particular,  we  are  interested  in  target  shapes  modelled 
by  homogeneous  polynomials  defined  on  elliptical  domains  or 
square  domains  with  rounded  corners 


2.2.  Conditional  symmetries  associated  with  this  model 


The  direct  method  (Clarkson  and  Kruskal  [5]) 
applied  to  a  second  order  PDE 


(9) 


consists  of  seeking  solutions  of  the  form 

E(x,y)  =  <&(x,y,F(z)),  z  =  z(x,y),  (x,y)e£l.  (10) 

•  z  =  z(x,y)  is  called  similarity  variable 

•  {( x,y )  :  z(x,y)  =  k}  are  called  similarity  curves 

-  after  substituting  (10)  into  (9),  we  require  that  the  result  to 
be  an  ODE  for  the  arbitrary  function  F  =  F(z) 

— ►  certain  conditions  are  imposed  upon  the  functions  <P,z  and 
their  partial  derivatives 

The  particular  case  [18] 


E(x,y)  =  F(z(x,y)) 


(11) 


consists  of  looking  for  solutions  depending  only  on  the  simi¬ 
larity  variable  z 

-  If  z  is  an  invariant  of  the  group  action  then  the  solutions  of 
the  form  (11)  are  as  well 

-  Assume  z  =  z(x,y)  is  s.t.  ||Vz||  ^0  on  2 


[5]  P.  A.  Clarkson  and  M.  Kruskal,  New  similarity  reductions  of  the  Boussi- 
nesq  equation,  J.  Math.  Phys.,  30  (1989),  pp.  2201-2213 


[18]  R.  Z.  Zhdanov,  On  conditional  symmetries  of  multidimensional  non¬ 
linear  equations  of  quantum  field  theory,  Symmetry  in  Nonlinear  Mathe¬ 
matical  Physics,  1  (1997),  pp.  53-61 


Assume  that  Young’s  modulus  takes  the  form 

E(x,y)  =  F(z(x,y))  (11) 

and  substitute  it  into  (8) 

— ►  we  get  the  following  relation 

F"(Z)  [Z^(WXX  +  VWyy)  +  2zXZy(l  ~  is)wxy  +  Zy(wyy  +  VWXX)\ 

~\~  F  (21)  \zxx(wxx  IS  yy')  -|-  2(1  is)  ZXyWXy  “|  \~Zyy(Wyy  “J-  IS  JXx) 

+22x(Aw)x  +  2zy(/\w)y\  +  F(z)(A2w)  =  1,  (12) 
which  must  be  an  ODE  for  F  =  F(z) 

—>  the  coefficients  of  the  partial  derivatives  of  F  are 
functions  of  z  only  (i.e.,  these  coefficients  are  also 
invariant  under  the  same  group  action) 

Denote  them  by 

Tl  (z)  =  Z2(wxx  +  VWyy)  +  2zXZy{\  ~  is)wxy  +  Z^  (_W  yy  +  VWXX), 

E  2^z)  -  ZXX(WXX  IS  yy')  ~\~  2(1  is)  ZXy  W  Xy  Zyy(  yy  “J-  IS  Uxx) 

+2zx(Aw)x  +  2zy(Aw)y , 
f~3(z)  =  A  2w. 

(13) 

— >  the  PDE  (8)  is  reduced  to  the  second  order 
linear  ODE 


r1{z)F”{z)  +  r  2(z)f’(z)  +  f3(z)f(z)  =  1  (14) 


2.2.1.  Data  and  parameter  invariant  under  the  same  group 


If  the  target  shape  is  invariant  under  the  same  group  action 
as  Young’s  modulus  then 

w(x,y)  =  G(z(x,y))  (15) 


-  substituting  (15)  into  the  relations  (13)  we  have 


fi  =  G"  (z2  +  z2) 2  +  G'  [(z2  +  vz2)zxx  +  2(1  -  v)zxzyzxy  +  (zy  +  VZ^Zyy]  , 
r2  =  2 Gr  (z2  +  z2)2  +  G"  {  [7z2  +  (y  +  2)zy]zxx  +  2(5  -  v)zxzyzxy 

+  [7Zy  +  (y  +  2)^]%y}  +  G'  {  (Az)2  +  2(1  -  v)Gxy  -  zxxzyy) 

+2  [zx(Az)x  +  zy(Az)y]}  , 


r3  =  G""{z 


X  +  )2  +  2  (3z 


X  +  Zy^XS  +  4^ 


zxy 


+  (^  +  3Zy)Z‘ 


yj^yy 


+  {3(Az)2  +  4(z2y  -  zxxZyy )  +  4  [zx(Az)x  +  z^A^);,] }  +  A2z. 

(16) 

-  the  coefficients  of  the  partial  derivatives  of  the  function  G 
must  depend  only  on  2,  i.e., 


where 

a2(z) 

ai(z) 

a2(z) 

a3(z) 

04(2) 

a5(z) 

clqG) 


|"i  =  a4G"  +  a\G\ 

r2  =  2  a4G"'  +  a2G"  +  a3G' , 

f3  =  a4G""  +  2a4G’///  +  cl5G"  +  a6G' , 


z2  +  z2, 

(z2  +  VZ2)ZXX  +  2(1  -  v)zxZyZxy  +  (z2  +  vzl)zyy , 


zxx 

,2 


+  2(5-1.)  ZXUyZXy  +  [7z2  +  (y  +  2)2 


yj  5 


7z2  +  G  +  2)z 

(Az)2  +  2(1  -  ^)(^xy  -  zxxZyy )  +  2  [zx(Az)x  +  zy(Az)y] 

(3^2  +  Zy^Zxx  “I-  ^ZXZyZXy  +  (^x  +  3zv)z 

3(Az)2  +  4(z 
A2z. 


'y  J^yyi 


2 

xy 


ZXX  Zyy )  +  4  [-Zx(^^.z)x  “I-  Zy(^^z')  y] 


yi  5 


-yy^ 


(17) 


From  the  first  relation  in  (17)  which  is  a  2D  eikonal  equation, 
we  get 

Z^ZxX  +  2  ZXZyZXy  +  Z^Zyy  =  «3  (  2  )  a'  (  2  )  , 

Zxx  =  a(z)a'(z)  -  zf  Zxy , 

Zyy  =  a(z)a'(z)  -  fzxy. 

The  last  two  equations  imply 

ZyZXX  -  2  ZXZyZXy  +  zlzyy  =  O'3  (z)  Oi  (z)  ~  OL^  (z)  .  (18) 

ZxZy 

Assume  that  there  is  a  function  /3  =  (3{z)  such  that 

Zxy  =  P(z)zxZy  (19) 

Indeed,  since  the  left  hand  side  in  (18)  depends  only  on  2, 
one  can  easily  check  that 

if  2  satisfies  both  the  2D  eikonal  equation  in  (17) 
and  (19),  then  all  the  functions  ai  —  a^z )  defined 
by  (17)  are  written  in  terms  of  a  and  f3 

— ►  the  problem  of  finding  z  is  reduced  to  that  of 
integrating  the  2D  eikonal  equation  and  the  PDE 
system 

Zxx  —  Cta'  -  (3z y 

<  Zxy  —  /3ZxZy  (20) 

,  zyy  =  aa>  -  Pz%. 

The  system  (20)  is  compatible  if  it  holds  the  following  relation 

aa"  +  a'2  —  3(3aa'  +  ol2  ( /3 2  —  (3'^j  =  0 


Denote  y  =  'a2.  In  this  case,  the  above  compatibility  condi¬ 
tion  can  be  written  as 

y"  -  3/V  +  2 y  (f32  -  (3')  =  0.  (21) 

On  the  other  hand,  if  the  function  [3  is  given  by 

af  ^  A"0) 

Kz)  =  -  w  (22) 

where  A  is  a  non-constant  function,  then  the  equation  (19) 
turns  into 

OO))*,,  =  0 


The  general  solution  of  this  equation  is  given  by 

\(z(x,  y ))  =  a(x)  +  b(y),  (23) 

with  a  and  b  being  arbitrary  functions.  Substituting  [3  from  (22) 
into  the  compatibility  condition  (21)  and  after  integrating 
once,  we  get 

//A'  +  2ii\"  =  fc,  (24) 

where  k  is  an  arbitrary  constant. 


Case  1.  If  k  7^=  0  then  after  integrating  (24)  and  substituting 
back  y  =  |a2,  we  get 


a. 


2 


2k\(z)  +  C  i 

A'2M 


(25) 


The  relation  (23)  implies  A \z)zx  =  a'(x),  and  A '(z)zy  =  6'(y) 


Substituting  these  relations,  (23)  and  (25)  into  the  2D  eikonal 
equation,  it  follows  that  a  =  a(x)  and  b  =  b(y)  are  solutions 
of  the  following  respective  ODEs 

a,2(x)  —  2  ka(x)  =  C2,  and  b'2(y )  —  2kb(y)  =  C3, 
with  C2  +  C3  =  Ci  (here  C*  are  real  constants) 


1 


The  above  ODEs  admit  the  non-constants  solutions 

a(i)  =  t[i2(I-C,)2-C2]  and  b(y)  =  ^[k2(y-C5)2-C3] 
and  so,  (23)  takes  the  form 


X(z(x,y))  =  ^  [O  -  C4)2  +  (y  -  C5)2]  -  E  (26) 
2  2k 


Notice  that  ^A  or  A  +  defines  the  same  function  3  as  the 

function  A  does.  Moreover,  since  the  PDE  (8)  is  invariant 
under  translations  in  the  (x,y)~ space,  we  can  consider 

A(^(x,  V ))  =  x2  +  y2  (27) 

If  VX  is  a  bijective  function  on  a  suitable  interval,  and  if  we 

denote  by  0  =  (Vx)  its  inverse  function,  then  z  written  in 
the  polar  coordinates  (r,6>)  (where  x  =  rcos(0),  y  =  rsin(0)) 
is  given  by 

z(x,y)  =  0(r)  (28) 

For  simplicity,  we  consider  0  =  Id,  and  from  that  we  get 

E  =  F(r),  w  =  (7(r),  z(x,y)  =  r  (29) 

Hence,  the  ODE  (14)  turns  into 


(g"  +  -G')  F"  +  ( 2G'"  +  ^AAq"  -  Ag')  F' 

V  r  J  y  r  r2  J 

+  ( G""  +  -G'"  -  Aq"  +  Ag')  F  =  1,  (30) 

\  ry»  E.  ry»  O  / 

which  can  be  reduced  to  the  first  order  ODE 

(G"  +  -G')  F'  +  (g"'  +  -G"  -  AgA  F  =  r  ~r°  +  (31) 

V  r  J  y  r  r2  J  2  r  r 

where  ro  e  [0, 1]  with  the  property  that 


7  = 


(rG" 


+  vG')  F'  +  (rG'"  +  G" 


F 


r—r0 


is  finite 


The  smoothness  condition  G"(0)  =  0  implies  that  the  equation 
(31)  can  be  written  as  [12] 

(g"  +  -G')  F'  +  ( G"'  +  —G"  -  4g')  F  =  -.  (32) 

V  r  /  y  r  r2-  J  2 

Case  2.  If  k  —  0,  similarly  we  get 

z(x,  y)  —  <t>(kix  +  k2y),  (33) 

where  k\  and  k2  are  real  constants  such  that  k\  +  k\  >  0 

In  this  case,  for  0  =  Id,  the  parameter  and  the  data  are 

written  as 

E  =  F(z)  w  =  G(z),  z(x,y)  =  k\x  +  k2y,  (34) 

and  the  ODE  (14)  turns  into 

G"(z)F"{z)  +  2  G"\z)F\z)  +  G""(z)F(z)  =  ,,2  *  2.2,  (35) 

(/cf  +  k^y 

with  {z\G"(z)  =  0}  the  associated  set  of  singularities. 

Integrating  the  above  ODE  on  the  set  {z\G"(z)  ^  0}  we  obtain 
that  Young’s  modulus  is  given  by 

(kix  +  k2y)2  +  Ci(k\x  +  k2y)  +  C2 
~  2 (k\  +  kl)2G"{klx  +  k2y )  ’ 

where  Ci  are  arbitrary  constants. 


[12]  E.  H.  Mansfield,  The  bending  and  stretching  of  plates ,  Cambridge 
University  Press,  1989 


2.2.2.  Data  and  parameter  invariant  under  different  groups 


Consider  two  functionally  independent  functions  on 
Q,  say  z  =  z(x,y)  and  v  =  v(x,y),  and  let 

w  =  H(v(x,y))  (36) 

be  the  target  shape 


-  In  this  case,  the  data  and  the  parameter  do  not 
share  the  same  invariance 

-  Similarly  to  the  above,  substituting  (36)  into  the 
relations  (13)  we  get 


r  i  = 


2 ZxZvVa 


+  ZyVxx^j 


r2  = 


r3  = 


=  [( ZXVX  +  ZyVy )2  +  V^ZyVx  ~  ZXVy )2 

+  z2vxx  +  2  ZX  ZyVXy  “I"  zyvyy  T  u  (  Zx^yy  ^-^x^yuxy  [  t'y'-'xx^ 

(v£  +  Vy)(zxVX  +  ZyVy)  +  V%ZXX  +  2  VXVyUXy  +  VyUyy 

-\-V  (VyUXX  -  2  VXVyZXy  +  V^Zyy)  +  2ZXVXVXX  +  2  (zxVy  +  ZyVX)vxy  +  2  ZyVyVyy 

~\~  ^Zxv X  H-  ZyVy)  (A^)]  H-  [-^XX^XX  +  2 ZXyVXy  “f"  ZyyVyy  IS  ( ZXXVyy  2 ZXyVXy 

H-  ZyyVxx  )  +  2x(A'y)x  +  Zy(Av)y] 

Cx  +  +  2  [(3vJ  +  +  4vxvyvxy  +  (>2  +  3^2  )v. 


j  uyy 


+  3^  +  Av'ly  +  3^y  +  2  vxxVyy  +  4^(A^)X  +  4vy(Av)y  +  A2v. 

(37) 


Recall  that  C’s  are  functions  of  2  =  z(x,y)  only.  Since  each 
right  hand  side  in  the  above  relations  contains  the  function 
H  =  H(v)  and  its  derivatives,  We  require  that  the  coefficients 
of  the  derivatives  of  H  to  be  functions  of  v.  It  follows  that  r 
must  be  constant  and  denote  them  by  7 \ 

Therefore,  the  last  condition  in  (37)  becomes 

A  2(w)  =  73,  (38) 


which  is  the  biharmonic  equation 


According  to  the  above  assumption,  we  seek  solutions  of  (38) 
that  are  functions  of  v  only 

Similarly,  we  get 

v(x,  y)  =  4/(r)  or  v(x,  y)  =  Ms(k±x  +  k^y)  (39) 
Thus,  for  T  =  Id,  the  target  shape  is  written  as 

w(x,y)  =  H(r)  or  w(x,  y)  =  H(k\x  +  k3y)-  (40) 

Since  2  and  v  are  functionally  independent,  we  get 


z(x,y)  =  fcix  +  v(x,y)  =  x2  +  y2  (41) 

or 

z(x,y)  =  +  y2,  v(x,y)  =  k±x  +  k2y  (42) 

One  can  prove  that  if  the  coefficients  7*  are  constant,  and  if  2 
and  v  are  given  by  (41)  or  (42),  respectively,  then  7!  =  72  =  0, 
and  73  7^=  0 

On  the  other  hand,  the  solutions  of  the  biharmonic  equation 
(38)  of  the  form  (40)  are  the  following 


ui(x,y)  =  -jr~zA  +  C\z2  +  C2  ln(z)  +  C322  ln(z)  +  Ca,  z=  J x2  + 

64  v 

and  respectively, 

73 


y‘ 


w(x,y )  = 


+  C3v  +  C4,  v  =  k^x+k^y 


24(/cJ  +  fc2)  2 

and  these  correspond  to  the  constant  Young’s  modulus 

.  N  1 
E(x,y)  =  — 

73 


(43) 


Notice  that  only  particular  solutions  of  the  biharmonic  equa¬ 
tion  have  been  found  in  this  case  (i.e.,  solutions  invariant  under 
rotations  and  translations) 

Since  this  PDE  is  also  invariant  under  scaling  transformations, 
which  act  not  only  on  the  space  of  the  independent  variables 
but  on  the  data  space  as  well,  it  is  obvious  to  extend  our  study 
and  to  seek  for  other  type  of  symmetry  reductions. 


2.3.  Equivalence  transformations  related  to  the  studied  model 

Consider  a  one-parameter  Lie  group  of  transformations  acting 
on  an  open  set  V  c  Q  x  W  x  S,  where 


W  is  the  space  of  the  data  functions  and 
S  is  the  space  of  the  parameter  functions 


given  by 

x *  =  x  +  sCO?  y? w-,  E)  +  0(s2) 

<  y*  =  y  +  £'n(x,y,w,E)  +  e>02) 

w*  =  w  +  y,  w,  E)  +  0(s2) 

K  E*  =  E  +  s^^x.y.w,  E)  +  0(e2), 

where  £  is  the  group  parameter.  Let 


(44) 


V  =  CO,  y,  w,  E)dx-\-r](x ,  y,  w,  i£)dy+</>0,  y,  w,  £;)<9w+^0,  2/ ,  w, 

(45) 


be  its  associated  general  infinitesimal  generator 


-  (44)  is  called  an  equivalence  transformation  associated  to 
the  PDE  (8)  if  this  leaves  the  equation  invariant 


-  the  equivalence  transformations  can  be  found  by  applying 
the  classical  Lie  method  to  (8),  with  E  and  w  both  considered 
as  unknown  functions 


following  this  method  we  obtain 


CO  ,y,w,E) 

T](x,y,w,E) 
c/)(x,y,w,E) 
,  ^(x,y,w,E) 


k±  +  k$x  —  k^y 
k2  +  k^x  +  key 

k3  +  kjx  +  key  +  (4 fc5  -  ko)w 
ksE 


where  ki  are  real  constants 


(46) 


-  the  vector  field  (45)  is  written  as  V  =  ^f=1  hVi,  where 

V!  =  dx,  V2  =  dy,  V3  =  dw, 

V4=-ydx  +  xdy,  V5  =  xdx  +  ydy  +  4wdw,  (47) 
V5  =  ydw ,  V7  =  xdw ,  Vs  =  -wdw  +  E8e • 


The  equivalence  transformations  related  to  the  PDE  (8)  are 
generated  by  the  infinitesimal  generators  (47) 

Thus,  the  equation  is  invariant  under  translations  in  the  x- 
space,  y-space,  w- space,  rotations  in  the  space  of  the  indepen¬ 
dent  variables  (x,y),  scaling  transformations  in  the  (x,y,w)~ 
space,  Galilean  transformations  in  the  (y,w)  and  (x,w)  spaces, 
and  scaling  transformations  in  the  (w,E)~ space,  respectively. 

-  Notice  that  the  conditional  symmetries  (§2.2)  represent  par¬ 
ticular  cases  of  the  equivalence  transformations 


-  Since  each  one-parameter  group  of  transformations  gener¬ 
ated  by  Vi  is  a  symmetry  group,  if  (w  =  G(x,  y),  E  =  F(x ,  y)) 
is  a  pair  of  known  solutions  of  (8),  so  are  the  following 


wEj  —  G(x  -  £i,y), 

w,(2  =  G(x,y  -  e 2 ), 

WG;  —  G(x,y)  +  £3, 

w(4)  =  G(x,y ), 

tt/5  =  e465^  (e_£5x,  e_£sy) 
w(6)  =  G(x,y)  +  £62/, 
w/7  =  G(:c,  y)  +  £7#, 

=  e~£8G(x,y), 


E™  =  F(x-euy), 

E^  =  F(x,y- e2), 

E&  =  F(x,y), 

EW  =  F(x,y), 

E^j  —  p  ^ e~£bx ,  e_£sy) 

E^=F(x,y), 

E™  =  F(x,y), 

E (8)  =  e£sF(x,y), 


where  7  =  xcos(£4)  +  ysin(e4),  y  =  -xsin(c4)  +  y  cos(e4),  and 
are  real  constants 


-  the  general  solution  of  (8)  constructed  from  a  known  one  is 
given  by 

u>(cc,  y )  =  eAes~£8G  (e~£5(x  —  k\),  e~£5(y  —  k2)^j  +  e4£5~£8£6y  +  e4£5-£8£7cc  +  e4£5_£s 
E(x,  y)  =  e£8F  (e-£5(5  -  fci),  e~£5(y  -  k2))  , 

where  k±  =  e\  cos(c4)+£2  sin (er4) ,  and  fc2  =  £1  sin(c4)-£2  cos(c4). 

-  The  equivalence  transformations  form  a  Lie  group  Q  with 
a  8-dimensional  associated  Lie  algebra  A.  Using  the  adjoint 
representation  of  Q,  one  can  find  the  optimal  system  of  one¬ 
dimensional  subalgebras  of  A  (more  details  can  be  found  in 
[3],  pp.  203-209).  This  optimal  system  is  spanned  by  the 
vector  fields  given  in  Table  1 

-  Denote  by  2,  I,  and  J  the  invariants  related  to  the  one- 
parameter  group  of  transformations  generated  by  each  vector 
field  Vi.  Here  F  and  G  are  arbitrary  functions,  (r, 6)  are  the 
polar  coordinates,  and  a,b,c  are  non-zero  constants. 


[3]  P.  J.  Olver,  Applications  of  Lie  Groups  to  Differential  Equations ,  Grad¬ 
uate  Texts  Math.  107,  Springer  Verlag,  New  York,  1986 


Generators 

Invariants 

£ 

II 

j? 

II 

ODE 

Vi 

z  =  y 

I  =  w 

J  =  E 

w  =  G{z) 

£  =  F(z) 

(50) 

Z  =  X 

I  =  w 

J  =  E 

w  =  G(z) 

£  =  F(z) 

(50) 

V4 

z  —  r 

I  =  w 

J  =  E 

iu  =  G(z) 

E  =  F(z) 

(30) 

V5 

z  =  y- 

X 

I  =  x~4w 

J  =  E 

w  =  .x4GY(z) 

£  =  F(z) 

(52) 

cV3  +  V4 

z  =  r 

I  =  w  —  cO 

J  =  E 

iu  =  c0  +  G(z) 

£  =  F(z) 

(30) 

Vs  +  cV8 

z  =  y- 

I  =  xc~4w 

J  =  x~cE 

w  =  x4_cG(z) 

£?  =  xcF{z ) 

(53) 

V4  +  cV8 

z  =  r 

I  =  ec0w 

J  =  e~c0E 

w  =  e~c0G(z) 

£  =  eceF(z) 

(54) 

V4  +  cVs 

— cO 

z  =  re 

I  =  r~4w 

J  =  E 

w  =  r4G(z) 

£  =  F(z) 

(56) 

V4  +  cX5 
+bVi 3 

— cO 

z  =  re 

T  -—4 

1  =  To  w 

J  =  r~cE 

w  =  r4_^G(2;) 

£7  =  rcF(z) 

(57) 

Vi  +  cV6 

z  —  y 

I  =  w  —  cxy 

J  =  E 

w  =  cxy  +  G(z) 

E  =  F(z ) 

(50) 

V2  +  cV7 

Z  —  X 

I  =  w  —  cxy 

J  =  E 

iu  =  cxy  +  G(z) 

£  =  F(z) 

(50) 

V 1  +  cVb 

z  =  y 

I  =  ecxit; 

J  =  e~cxE 

w  =  e~cxG{z ) 

£?  =  ecxF(z ) 

(58) 

V2  +  cVb 

z  =  X 

I  =  ecyiu 

J  =  e~cyE 

iu  =  e~cyG{z ) 

E  =  ecyF(z ) 

(58) 

Table  1 


— ►  To  reduce  the  order  of  the  PDE  (8)  one  can  integrate  the 
first  order  PDE  system 


f  C (x,y,w,E)wx  +  ri(x,y,w,E)wy  =  c/)(x,y,w,E) 

I  C (x,y,w,E)Ex  +  rj(x,y,w,E)Ey  =  ij)(x,y,w,E), 

that  defines  the  characteristics  of  the  vector  field  (45) 


(49) 


The  reduced  ODEs  listed  in  Table  1: 


•  The  invariance  of  the  equation  (8)  under  the  one-parameter 
groups  of  transformations  generated  by  Vi,  V2,  V\  +  cV6  and 
V2  +  CV7,  respectively,  leads  us  to  the  same  ODE 


F\z)G\z)  +  2  F\z)G"\z)  +  F(z)G""(z)  =  1,  (50) 

with  the  general  solution 


F(z) 

on  the  set  { z\G"(z )  ^  0} 


z2  +  ciz  +  c2 

2  G"(z) 


(51) 


•  The  invariance  under  the  scaling  transformation  generated 
by  the  vector  field  V5  yields  the  reduced  ODE 


G"  (z2  +  l)2  -  6 z(z2  +  1  )G'  +  12 (z2  +  v)G 


F' 


+  2 


(z2  +  l)2  G'"  -  5z(z2  +  1  )G"  +  3(4z2  +  v  +  1  )G'  -  12 zG 


F' 


+ 


(z2  +  l)2  G""  -  4 z(z2  +  1  )G'"  +  4(3^2  +  1)G"  -  24 zG'  +  24G 


F  =  1. 
(52) 


•  The  following  ODE 


(z2  +  l)  G"  +  2(c  —  3)z(z2  +  1)G'  +  (c  —  3)(c  —  4)(z2  +  v)G  F" 

+  1 2  [z2  +  1  )2  Gm  +  2(2c  —  5)z(z2  +  1  )G"  +  2(c  -  3) [z2(c  -  4)  +  v(c  -  1)  -  1  ]G' 

— 2(c  —  3)(c  —  4-)zG}  F1  +  j(z2  +  l)2  G"n  +  2(c-  2)z(z2  +  1)G///  +  [(c  -  3)(c-  4) 

— 2(c  -  2)  +  z/c(c  -  1)]  G"  -  2(c  -  4 )(c  -  3)zG'  +  2(c  -  4)(c  -  3)G}  F  =1 

(53) 


is  obtained  in  the  case  6.  The  reduced  equation 


,,  V  ,  VC 

G  +  —G  +  —  G 


F"  + 


2  cr  +  V-±^G"  +  _  cFl+F}G 

rp  p2  p3 


F' 

(54) 


+ 


G""  +  —G'"  + 


c2v  —  1  1  —  c2(2v  +  1)  .  2c2  (v  +  1) 


-G  + 


G  + 


-G 


F  =  1 


is  related  to  the  case  7.  This  can  be  written  as  the  first  order 
ODE 


.2  „2 


„  v  .  vc \  .  (  ...  1  „  c2v  —  1  ,  c2(  1  T  v)  \  r"  —  rx  ■-f 

G  +  -G  +  —  G  )  F'+  Gm  +  —G  -\ - - — G - v  J  J  G  )  F=  0  1  1 


r - 


2  r  r 

(55) 


where  ro  e  [0, 1]  with  the  property  that 


* 

7  = 


F'  (  rG"  +  vG'  +  -G)  +  F  (  rG'"  +  G"  +  — - — G'  -  C  ^  j~^G 


|r=ro 


is  finite 


•  In  the  cases  8  and  9,  after  the  change  of  the  variable 
2  =  exp(t),  the  reduced  ODEs  are  the  following 


{(c2  +  1  )20'  +  (c2  +  1)0  +  7)0  +  4  [03 c2  +  1)  +  c2  +  3]  G)  F" 

+  {2  (c2  +  1)20"  +  (c2  +  1)0  +  19)0'  +  2[16  +  (c2  +  1)(3  z/  +  13)]0 
+  80  +  7  )G}F'  +  {(c2  +  1)20'"  +  12(c2  +  1)0" 

+  4(5c2  +  13)0'  +  960  +  64 G}F  =  1, 


(56) 


and  respectively, 


|  (c2  +  l)2  +  ^ - 1-  [cO  +  7)  —  26](  +  ^ - ^  [c3 ( 1  +  3z/)  —  c2z/& 

+c(i/  +  3)  -  &]G}  F"  +  |2(c2  +  1)20"  +  (i  +  c)  [c(i/  +  19)  -  4 6] O' 


+  2 


— —  -(-  z/'fr2  -(-  c2(3 v  13)  —  Abc{y  -(-1)  —  12 — f-  3 v  29 


G' 


+  P-4 


)  [2c(v  +  7)  +  6(1/  -  5)]g|  F'  +  {(c2  +  1)2G 
+2  (c  +  i)  (6c  -  6)G"'  + 

+  (®  -  4)  [16c  +  6(i/  -  5)]G'  +  2  (4  -  0  I8c  +  K"  -  3)]g|  F  =  1. 


4  +  _(„  _  17)  +  20c2  -  6c (c  +  7)  +  J/62  +  52 

c2  c 


O 


(57) 


•  In  the  cases  12  and  13  we  get  the  same  equation 

(O'  +  uc2G )  F"  +  2  (O"  +  z/c20)  F'  +  (O'"  +  z/c20')  F  =  1,  (58) 

with  the  general  solution  given  by 

z2  +  Ciz  +  C2 


F(z)  = 


(59) 


G"(z)  +  uc2G(z ) 

on  the  set  {z\G"(z)  +  uc2G(z)  ^  0},  where  C\  and  C2  are 
arbitrary  real  constants 


2.4.  Examples 


— ►  Since  the  data  w  is  the  function  that  models  the  target 
shape  of  a  car  windscreen,  we  seek  a  data  with  relevant  phys¬ 
ical  and  geometrical  properties,  such  as  smoothness  and  a 
positive  curvature  graph  at  least  in  the  center  of  Q,  for  which 
the  boundary  condition 

w(x,y)\an  =  0  (6) 

is  satisfied 


Moreover,  if  there  is  no  free  rotation  of  the  plate  around  the 
tangent  to  <9Q  then  the  condition 


is  required. 


dw 

dn 


3Q 


o  (7) 


Remarks.  Assume  that  =  {(x,y) \  z(x,y)  =  k}  is  the  k- 
level  set  of  the  function  z  (here  k  being  a  non-zero  constant) 
and  ||Vz||  >  0  on  Q.  If  the  target  shape  is  given  by 


w(x,y)  =  a(x,y)G(z(x,y )) 

where  a  =  a(x,y)  is  a  suitable  function  according  to  Table  1, 
then  the  boundary  conditions  (6)  and  (7)  are  equivalent  to 


G(k)  =  0  and  G\k)  =  0 

Therefore,  the  data  might  have  the  form 

w(x,y)  =  a(x,y)  ( z(x,y )  -  k)  H(z(x,y)) 

This  corresponds  to  the  case  when  the  data  and  Q  are  invari¬ 
ant  under  the  same  symmetry  reduction 

In  our  case,  this  can  be  applied  to  rotational  symmetries 


— ►  For  scaling  invariance,  we  have  to  incorporate  the  non¬ 
invariant  boundary  conditions  in  invariant  solutions 

-  the  problem  on  elliptical  domains  and  on  square  domains 
with  rounded  corners 

For  instance,  the  class  of  target  shapes  of  the  form 


w(x,y )  =  zm(x,y)  -  km 

where  m  >  la  natural  number,  satisfies  (6) 
In  this  case,  we  have 


dw 

dn 


on 


=  mkrn~1\\'Vz\ 


dQ 


If  this  quantity  is  small  then  the  condition  (7)  is  almost  sat¬ 
isfied  (i.e.,  there  is  a  small  rotation  of  the  plate  around  the 
tangent  to  <9Q) 

In  the  following  examples,  we  assume  that  the  glass  Poisson 
ration  v  =  0.5 


E.g.  1.  Rotational  invariant  data  and  parameter 
Consider  the  target  shape  of  the  form  [14] 

w(x,  y)  =  G(r)  =  (r  -  l)2  (2 r  +  1)  ,  r  =  \/ x2  +  y2, 

o  v 

defined  on  the  unit  disc  which  satisfies  the  boundary  conditions 
(6)  and  (7).  Since  G'( 0)  =  0,  the  reduced  ODE  is  (32)  and 
this  has  a  singularity  at  r  =  |.  Since  E  >  0,  we  consider  the 
constant  of  integration  Ci  =  1,  and  so, 

E(x,  y)  =  F(r)  =  ^  +  (5r  -  3)'^ 

The  PDE  (8)  is  elliptic  for  r  e  (|,  |),  hyperbolic  for  r  e 
[0,  |)  |J  (|,  l] ,  and  parabolic  if  r  =  |  or  r  =  respectively. 
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Fig.  1 


[14]  D.  Temple,  An  inverse  system:  An  analysis  arising  from  Windscreen 
manufacture,  M.  Sc.  Thesis,  Department  of  Mathematics,  Oxford  Univer¬ 
sity,  Oxford,  United  Kingdom,  2002 


E.g.  2.  Particular  target  shapes  on  rounded  square  domains 

a)  Suppose  the  Lame  oval  <9Q  =  {(x,y) \  x2n  +  y2n  =  1}  is  the 
boundary  of  the  domain  (here  n  >  2  is  a  natural  number).  For 

w(x ,  y)  =  ( x2n  +  y2n)m  -  1,  (60) 

m  >  1  being  a  natural  number,  the  equation  (8)  is  elliptic 
on  Q  — {(0,0)}  and  parabolic  in  (0,0) 

These  target  shapes  are  invariant  under  V5  +  cVq  +  (4  -  c)V 3, 
where  c  =  4  -  2 ran.  For  1  >  0  or  1  <  0,  the  functions  (60) 
can  be  written  as 

w(x,y)  =  x2mnG(z)  -  1,  G(z)  =  (l  +  z2n)m ,  z  = 

X 

According  to  the  case  6  in  Table  1,  the  associated  Young’s 
modulus  has  the  form 

E(x,  y )  =  x4~2mnF(z),  z  =  -. 

X 

Since  w(x,  y)  =  w(y,  x)  =  w(—x,  y)  =  w(x,  —y)  =  w(—x,  —y) , 

Young’s  modulus  also  shares  these  discrete  symmetries.  Thus, 
the  reduced  ODE  (53)  can  be  integrated  for  2  e  [0,1].  For 
n  —  2  and  m  =  1,  the  data  is  a  solution  of  the  biharmonic 
equation  and  £;  =  48_1. 

For  n  —  3  and  m  —  1,  the  data  and  the  numerical  solution  F 
satisfying  F( 0)  =  0.002  and  F'( 0)  =  0  are  given  in  Fig.  2. a. 


Fig.  2.  a 


b)  Assume  that  <9Q  =  {(x,y)\  x2n  +  y2  =  1},  where  n  >  1  is  a 
natural  number.  Consider  the  following  class  of  target  shapes 

w(x,y)  =  x2n  +  y2  -  1, 

is  invariant  under  the  vector  field  V2-\-2Vq.  Hence,  the  equation 
(8)  is  reduced  to  the  ODE  (50) 


For  n  —  3  and  m  =  1,  the  associated  Young’s  modulus  is 


E(x,y)  =  F(x)  = 


x2  +  C  ix  +  C2 

2(30x4  +  l) 


Since  E  >  0,  we  can  set  Ci  =  0  and  C2  =  1  (see  Fig.  2.b). 
The  PDE  (8)  is  elliptic  on  Q.  —  Oy  and  parabolic  on  the  y-axis. 


Fig.  2.b 


E.g.  3.  Particular  target  shapes  on  elliptic  domains 
Consider  the  data 

on  the  elliptic  domain  Q  =  {(x,y) |  ^  ^  <  1},  where  m  >  1 

is  a  natural  number.  These  target  shapes  are  invariant  under 
Vs  +  cV s  +  (4  -  c)  V3 ,  where  c  =  4  -  2m. 

The  PDE  (8)  is  elliptic  on  Q  — {(0,0)}  and  parabolic  in  the 
origin.  For  x  >  0  or  x  <  0,  (61)  can  be  written  as 

w(x,y)  =  x2mG(z)  -  1,  G(;z)  =  (  —  +  ^  )  ,  z  — 

\az  bz  J  x 

In  this  case,  we  look  for  solutions  to  (8)  of  the  form 

E(x,y )  =  x4-2mF(z),  z  =  — 

x 

If  m  =  2,  the  data  is  a  solution  of  the  biharmonic  equation,  and 
the  related  Young’s  modulus  is  E  =  24  (a-4  +  6~4)  +  16a~26-2. 

If  m  >  3,  the  reduced  ODE  is  (52).  For  m  =  3,  the  data  and 
the  numerical  solution  F  of  (52)  satisfying  F( 0)  =  0.001  and 
F'( 0)  =  0  are  plotted  in  Fig.  3 


Fig.  3. 
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